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Abstract 

We use the chain of simple heuristic expedients to obtain perturbative and ex- 
actly solvable relativistic spectra for a family of two-fermionic bound systems with 
Coulomb-like interaction. In the case of electromagnetic interaction the spectrum 
coincides up to the second order in a coupling constant with that following from 
the quantum electrodynamics. Discrepancy occurs only for S-states which is the 
well-known difficulty in the bound-state problem. The confinement interaction is 
considered too. 

PACS number(s): 03.65.Pm, 03.65. Ge, 12.39.Pn 

1 Introduction 

A nonrelativistic two-body problem reduces completely to the single-body problem with 
a central potential. The single-body Schrodinger equation becomes the two-body one (for 
the centre-of-mass frame of reference) if to replace the particle mass by the reduced mass, 
and to understand the radius-vector r as the relative position vector. 

In the relativistic case the relation between the single- and two-particle problems is not 
so transparent. There are few reasons for that. First, a spin appears on a physical stage. 
The existence of spin diversifies properties of both interacting particles of matter and fields 
mediating this interaction. Thus even in the single-particle problem various relativistic 
wave equations such as the Klein-Gordon, Dirac, Duffin-Kemmer and other equations 
exist, and they can involve scalar, vector, tensor and other interaction potentials. Second, 
there exist different approaches to the relativistic two- and few- body problem. The 
most profound approaches are based on the quantum field theory (QFT), especially on 
the perturbative QFT. They lead to complicated integral equations (such as the Bethe- 
Salpeter equation [0, 0), coupled sets of differential equations (such as the Breit equation 
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H Q]), or their higher order differential reductions (such as the Fermi- Pauli equation 
0). Other approaches such as quasipotential 0-0 or ones based on the relativistic 
direct interaction theory (RDIT) f|-0 are semi- or purely phenomeno logical. They 
manifest a general structure of relativistic potentials and wave equations while features of 
concrete interactions must be brought from other sources (for example, from the classical 
or quantum field theory). 

I. Todorov has observed a simple way how the Klein-Gordon equation can be trans- 
formed into the quasipotential equation describing the scalar and vector interaction of two 
spinless particles 0,0. In the case of Coulomb-like interaction this equation gives the 
spectrum which agrees with QFT result up to a 4 -terms of coupling constant expansion. 
The generalization for an arbitrary field-type interaction (including higher rank tensor 
interactions) arises naturally from the Schwinger source theory and the Fokker for- 
malism ||, and results in some RDIT models [TT]-[T3|]. In the present paper we construct 



relativistic wave equations appropriate for the description of the field-type interaction of 
two particles with spin \. 

At the beginning we summarize results concerning spinless system. Namely, in Section 
2 we formulate the family of relativistic wave equations which describe the scalar, vector 
and gravitational (i.e., tensor) interaction of two scalar particles. These equations have a 
common effective single-particle structure. In the case of Coulomb-like interaction they 
are exactly solvable. The mass spectra coincide with that obtained algebraically in Ref. 
13| , and agree up to a 4 with known QFT results. The only disagreement exists for 



S-states. 

Then we modify the wave equations in order to describe two-fermion systems. The 
spin-orbital, spin-spin and tensor corrections to the scalar, vector and gravitational in- 
teractions are brought from QFT§ (Section 3). In Section 4 by means of an appropriate 
rescaling of r spin interaction is included into the eigenstate problem as small correction to 
the effective Coulomb Hamiltonian. Calculations with the first order perturbation theory 



(Section 5) reproduce the QED muonium spectrum up to a 4 [15| and give a generalization 
for the cases of scalar and gravitational interactions. 

Spin potential corrections depend on r = |r| as 1/r 3 . Due to properties of matrix 
elements with Coulomb bound states we transform spin corrections into 1/r 2 terms in 
such a way that perturbative spectrum (up to a 4 ) does not change. In this form the 
problem appears exactly solvable, which is shown explicitly in Section 6. 

The Todorov construction was proposed in the case of Coulomb-like interactions. Here 
(in Section 7) we approve this recipe for a system with confining interaction, and make 
an appropriate modification to account spin effects. 



2 Spectra of systems of two spinless particles 

Let us consider the stationary Klein-Gordon equation for particle of the rest mass m in 
the scalar potential V s (r) and the vector one V v (r) depending on r = |r|: 

Att(r) + ([E - K(r)] 2 - [m + V s {r)f) *(r) = 0. (1) 

1 We note that two-particle spin-orbital interaction can be recovered completely from the single-particle 
one [flif but this is not concerned with spin-spin and tensor terms. 
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Following the Todorov's observation in the quasipotential approach || [/] one can construct 
the appropriate two-particle wave equation by the following substitution: 



M 2 — m? — ml mim 2 
E — > Em = 2M ' m ^ m ^ = ^-' ( 2 ) 

where m a is the rest mass of ath particle, and M is the total mass of the system, i.e., 
the energy in the centre-of-mass frame of reference. One obtains the Schrodinger-like 
equation: 

A*(r) + [Q M -[/ M (r)]tf(r) = 0, (3) 

where 

Qm = E 2 M - m 2 M = J^[M 2 - m 2 + ][M 2 - m 2 _], (4) 
is the on-shell value of the relative momentum squared as a function of M [[7], 

U M (r) = 2[m M V s (r) + E M V v (r)] + V 2 (r) - V 2 (r) (5) 

is the effective potential, and m± = mi ± m 2 . In the nonrelativistic limit the equation 
(|3]) becomes the usual Schrodinger equation with the nonrelativistic potential V(r) = 
V B (r) + Vy(r) where r is the distance between particles. 

The effective potential Um{j") depends of the total mass M. Thus the equation (||) is of 
the quasipotential type, and problems can occur with the consistent quantum-mechanical 
treatment |7|]. 

In the case of Coulomb-like interaction 

V a (r) = -a B /r, (6) 
K(r) = -a v /r (7) 

with the coupling constants a s > and a v > the equation (H) is exactly solvable. 
Moreover, for the case of electromagnetic (vector) interaction the corresponding spectra 
coincide (except for the ground S-states) with those following from the scalar quantum 
electrodynamics in the second order of a coupling constant |6|, |5|] . 

The simple Todorov recipe embraces the cases of scalar and vector interactions (and 
their superposition). The generalization to the cases of higher-rank tensor interactions can 
be built on the base of the family of RDIT models (known as time-asymmetric) ]TT], [12 . 
This leads to the Schrodinger-like (quasipotential) equation (|3[) with 

U M (r) = -2m M f(X)- + ^, (8) 

where A = Em/ttim, and the function /(A) (such that /(l) = 1) as well as the constant 
7 depend on the tensor nature of interaction |12], |IBf . For example, for the scalar, vector 
and gravitational interaction (or another second-rank tensor interaction, for example, the 
strong gravitation) we have: 

/b(A) = 1, 7s = 1, (9) 

/v(A) = A, 7v = -l, (10) 

/ g (A) = 2A 2 -1, 7g = -6. (11) 
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The superposition of these interactions can be considered by means of the superposition of 
the functions and constants (|9|)-(|TT|). Also we note that for scalar and vector interactions 
(and their superposition) the effective potentials (§) (from the time-asymmetric models) 
is identical to (|5|) with (||) or/and (|7|) (from the Todorov recipe). 

The mass spectrum following from the equation (^|), @, (||) can be presented in the 
implicit form: 

M 2 = m\ + ml + 2m 1 m 2 X, (12) 



where A is a positive solution of the equation: 

1 - A 2 a 2 



/2(A) ^' 

and 



(13) 



v = n r + \ + v /(£+i) 2 + « 2 7 (14) 

is effective "principal quantum number"; here n r = 0, 1, ... is the radial quantum number 
and i = 0, 1, ... is the angular quantum number. 

For the scalar, vector and gravitational interactions the equation is solvable: 



A s = yf\ - a 2 /u 2 , (15) 



A v = 1 Jl + a 2 /u 2 , (16) 



1 / . V 2 V 



g 2 V / 2V a 2 a 
Approximately, with accuracy up to a 4 , we have 



4- — + -J8 + u 2 /a 2 . (17) 



m r a 2 m r a A 



M m m + — + 



2n 2 2n 4 



4 4m, 



m r a A 7 
+ ' 



n 



3 21+V 



where /' = df/dX, m r = m 1 m 2 /m + is the reduced mass, and n = n r +£+l is the principal 
quantum number. 



3 Two-fermion systems: including spin effects 

The weakly relativistic system of two fermions interacting via scalar or/and vector field 



can be described by the generalized Breit-Fermi Hamiltonian |Tf6| . |T7j| . Besides the nonrel- 
ativistic Coulomb Hamiltonian, it includes relativistic kinematic terms, spin-independent 
and spin-dependent corrections to the interaction. Some of these terms are singular and 
can be taken into account as perturbations only. 

Here we do not consider the Breit-Fermi Hamiltonian. Instead, we modify the Todorov 
recipe in order to describe the spin effects in two-fermionic systems. For this purpose we 
need only a spin-dependent part W of the Breit-Fermi Hamiltonian. For the scalar and 
vector interaction it is |T6|, \Fi 
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+ 



2 

mf 



+ 



m\ 



+ 



m\m2 



(T'2 



CTi • <X 2 AKW 



(20) 



where L = — i r x V is the orbital momentum operator, a a is the spin operator acting on 
the ath particle spin variable, and T = 3(<Ti • n)(<x 2 • n) — <T\ • er 2 is the tensor operator. 
In the case of gravitational interaction (with the nonrelativistic potential V g (r) = —a g /r, 
where a g = Gmim2, and G is the gravitational constant) we have fllSj : 




m\ 



+ 



mim 2 



c:v„ 



c 2 



4mim 2 r 



3mim 2 



a 1 ■ <r 2 5(r). 



(21) 



Now in order to construct a two-fermion equation we replace a nonrelativistic poten- 
tial V by V = V + W in the effective potential U M (Eq. (g) or (g)). The resulting 
quasipotential equation is not solvable, and we should apply some approximate method. 

In the case of Coulomb-like interaction the spin term W is meant to be small as to 
compare to the nonrelativistic potential V(r) = —a/r. Thus we can modify approximately 
the effective potential M) as follows: 



U M {r) - Um{t) « U M {r) + 2m M f(X)W. 



(22) 



Now one can account the spin correction by means of the perturbation method. In 
so doing we note the following. First, the original (non-perturbed) equation (|3]) fails to 
describe correctly S-states. Thus we will neglect 5-functional terms in W (i.e., the last 
term in r.h.s. of Eq. (f2~l|), and the last term in r.h.s. of Eq. (|20| ) in the case (0)) as they 
contribute in S-states only. Second, the modified equation (as well as the unperturbed 
one (Q)) is quasipotential but not the true Schrodinger equation. Thus it needs some 
minor reformulation to be tractable within the perturbation method. 



4 Spin corrections to the Coulomb-like interaction 

In the case of Coulomb-like interaction the two-fermion wave equation reads: 



p 2 - 2m M f(X) \--W 



a 



7- 



a 



Q 



V = 0, 



where p = — i V. Using the substitution 

r = p/R M , p = R M 7r with R M = am M f(\) 



we present the equation (|23|) in the dimensionless Hamiltonian form: 

HV = etf. 

Here 

H = ff(°> + 



(23) 

(24) 

(25) 
(26) 
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is the total Hamiltonian, 

H (0) = 1^2 _ I (27) 

2 p 

is the basic Coulomb Hamiltonian, 

# (1) = a 2 r, r = -L + is(n) (28) 

2p z p 6 

is a perturbation since a is considered as a small parameter, and 

Qm A 2 - 1 



2R 2 M 2a 2 f 2 (X) 



(29) 



is a dimensionless energy (i.e., a spectral parameter). 

With sufficient accuracy (i.e., up to terms ~ 0(a)) the last term of T is equal to 
W I [o?Rm) but does not depend on M. This is provided by use of approximate equality 
M m m + in small terms. The general form of the operator £ acting on angular and spin 
variables is: 

£ = [(£ - S 2 )L ■ a + + 2r]SL ■ <x_ + ((1 - S 2 )T] /16, (30) 

where a± = <7\ ± cr 2 , 5 = m_/m + , and the constants £, 77 and £ for scalar, vector and 
gravitational interactions are defined as follows: 

(31) 
(32) 
(33) 



& 


= -1, 


= 1, 


Cs = o, 




= 3, 


Vv = -1, 


Cv = l 




= 7, 


= -3, 


Cv = l 



5 Basic states and first-order perturbation theory 



The basic Hamiltonian (|27 ) commutes with operators of orbital angular momentum L, 
total spin S = |o+, total angular momentum J = L + \cr + and parity P. In order 
to write down the basic eigenf unctions \E f ^°- ) (p) we use the angular "bispinor harmonics" 
if l (n) (i = A, 0, — , +). In 2x2 matrix representation they are |H]]: 



-1 

1 



(34) 



2j(j+i; 



2(j + 1)(2j+3) 



2j(2j-l) 



(j-p+l)( J+ p)Yj 



v/(j+//+l)C7'-At)^ 



(j-^+l)(j-^+2)^ + - 1 1 -J( J+/i+ i)( J - /i+ i)Y£ 1 



(35) 



(j+/i+l)(j-/i+l)^ 



3+1 



(j +/i+ l)(j+/i+2)^ +;i 



■3+1 



(j+Zi-^O'+Zi)^ 



/x-l 



(36) 
(37) 
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where ^(n) (// = —£,...,£) are the spherical harmonics depending on the direction 
n = r/r. The bispinor harmonics form an orthonormal set, in the sense that (i\k) = 

J dnTr(<£>J = 8^ where the integrations are taken over the entire solid angle. 

The bispinor harmonics posses the following properties (besides those due to properties 
of the spherical harmonics): 



L-*<p A = + 1)/, L atp = -(j + 2)<p-, 



where the components of the vector operator a are the Pauli matrices. 

The action of spin operators on the bispinor harmonics is as follows: cr\^> = (Ttp, 
cr 2 {p = {pcr T . We note that tp A is antisymmetric and (p 0,± are symmetric matrices. Then 
it follows from this and Eqs. (ESI) that ip A and <y9°' ± satisfy the following equalities: 



JV = j'O' + iM 3 = o,i,..., 

J 3 (p = fj,<p, fj, = 

LV = £(£+l) v \ £ = 



sy = s(s + i)(p\ s = 

p<p A fi = (-)V' , P<p T = R j+ V T , 

±L-<T + (p A = 0, ±L-<r + (p- = -(j+2)<f- 

iL-o- + ^° = -<p°, \-L-<T + ^ + = (j - 

iL-<r_V° = JjU + 1V M ; \^-<T-^ = 0, 

cr l ■ cr 2 (p A = -3(p A , er l ■ er 2 ip 0,T = y? 0,T , 



3, 


i = A,0, 




i = T, 


0, i 


= A, 


1, i 


= 0,=F, 



(39) 



lm A n It - o V + 7+2 - 



Now one can choose four independent basic eigenfunctions (p) (i = A, 0, — , +) of 
H^> as follows: 

= -u n Mv A '\nl $?W = ^+i(p)^(n), (40) 

where u n ^(p) is a solution of the radial Coulomb problem 

H t u n , t (p) = e^u n j{p) (41) 

with the effective Hamiltonian 

and the dimensionless eigenenergy 

= -l/(2n 2 ), n = l,2,... (43) 
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We note that the basic eigenfunctions ^ao(p) have the parity P = (— and fy^\p) 
have the parity P = (— ) J ' +1 . The function ^/^(p) describes the singlet (s = 0, £ = j) 
state while ^^(p) correspond to triplet (s = 1, i = j,j ± 1) states. 

Let us calculate the first-order correction e^ 1 ' to the dimensionless energy e ~ + 
a 2 e^\ The total Hamiltonian H = + a 2 T commutes with operators of parity P and 
total angular momentum J = L + \cr+. One can choose the wave functions \I/(p) as the 
eigenfunctions of J 2 , J3 and P. Thus they can be spanned onto states \E^q if P = (— 

or onto if P = (— ) J=Fl . In the each parity case zero-order eigenvalues are twice 
degenerated. Thus in the first order of perturbation theory we have 



.(1) 
■■(i,k) 



r« + r fcfc ± v (r^ — r 



kk) 



+ 4F: 



(44) 



with 



A, k = if P = and i 



k = + if P 



T = [Tik] is defined as follows: 



r = [r ik ] = [(i\r\k)] 

Taking (|28|) and ([E]) into account we have: 



J drTi ^*J(r)r* fc (r; 



) J , where the matrix 



(45) 



if p 



V, and 



2 
2 



(i+i|p- 2 b'+i) 



(A\V\A) (A\m 
(0\E\A) (0|E|0) 







(46) 



if P = (-)^' +1 , where 



and 



In particular, 



0'-i|p- 2 |j-i> 

(j+l\p^\j+l)(A\E\A) (j+l|p- 3 |j-l)(A|S|0) 
(j-l\p-3\j+l){0\Z\A) (j-l|p- 3 |j-l)(0|S|0) 

(i\E\k) = J dnTi (<pt(n)E<^(n)) 
(f|p s |£) = / dpu n/l (p) p s u n>e (p). 



1 



\P~ 



Ir 1 ' n 3 (£ + i)' '"' £(£+1) 

(£+l\p- 2 \£-l) = 0, (£+l|p~ 3 |£-l) = 0. 



(47) 

(48) 
(49) 

(50) 
(51) 



The relations (|50[) are well known in literature (see J20J or and fl51|) can be calculated 
by means of formulae given in p], chap. Mathematical Supplements, § f]. 



8 



Using (|50D , fl5ll) and calculating the matrix elements (z|£|fc) by means of Eqs. 
fl33|), (|39|) one obtains the matrix T and then the corrections to the dimensionless 
energy. Then, using fll~2|), (f29|) and expanding the total mass M in a one obtains the 
first-order mass spectra (i.e., with accuracy up to a 4 ). 

Due to the relations (|50| ) and (|5l| ) the matrix T is not diagonal if P = (— ) 3 . Thus the 
correspondent first-order states are the mixture of singlet (s = 0, £ = j) and triplet (s — 1, 
£ = j) states. In the P = (— case T is diagonal, and the triplet (s = 1, £ = j ± 1) 
states does not mix. Thus it is convenient to classify the first-order mass spectra by j 
and £. These spectra can be obtained from Eq. (0) by the following substitution: 



7 ->lf + <l>(t,j) 



(52) 



where the function <p(£,j) depends on both a spin state of the system and the tensor rank 
of mediating field. We have: 



8£{£ + 1) 
1 + 5 2 



1 + 5 2 ± Ml + 5 2 ) 2 + m 2 £{£ + 1) 



1 + 5 2 



4(£+iy 

1 



3, 

3 + 1, 
3-1, 



(53) 



1 

1 



+ 1J 

1-5 2 
2(2£-l)' 
1-6 2 



l±Jl+45 2 £(£+r 



3 + 1, 



(54) 



■1) 
3 



3 

27 



2(2£ + 3)' 
8 2 



A5 2 



+ r 



2(2^-1)' 
1-S 2 



+ 



3 + 1, 



3-1- 



(55) 



_ 2{£ + l) ' 2(2£ + 3)' 

The Eqs. (|18]), (|52"D and (|54"1) reproduce the muonium spectrum and (if mi 
m2 = Tn) the positronium spectrum M. 



6 Solvable simulation of first-order mass spectra 



Solving the Schrodinger equation fl2"5] ) perturbatively is due to the fact that spin inter- 
action term in the operator ([28|) depends on p as p~ 3 . Below we construct some exactly 
solvable model which reproduces the spectrum of perturbation theory. 
Let us modify the operator (|28| ) as follows 

r_f^, (56, 
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where 

Z(n) = 7 + 2{£(n)/L 2 } ordered . (57) 

The operator Z acts on angle and spin variables. It is not defined on states which contain 
the S-wave, but we refuse these states from the very beginning. On other states Z is 
supposed to be Hermitian. Thus it must be somehow ordered if E(n) and L 2 do not 
commute. 

It is easy to examine by means of Eqs. Q5PJ)-(|5T|) that in the first order of pertur- 
bation theory the Hamiltonian H = + a 2 G has the same spectrum as the original 
Hamiltonian H. This result does not depend on the ordering rule used in Z. 

Below we show that the new Schrodinger equation is exactly solvable. Of course, 
the exact solution and corresponding spectrum depend on the ordering rule. One can 
consider, for example, the following rules: 

{£/L 2 } ordered = i(S|L|- 2 + |L|- 2 S), (58) 
{S/L 2 } ordered = |L| 1 S|L| 1 (59) 

oo 

{S/L 2 } ordercd = |dte^ L2 Se-5 L2 , (60) 
o 

where |L| = \f\?. The last rule is inspired by the Feynman representation of an inverse 
operator: A~ x = J °° dt exp(-tA) . 

The radial reduction of the Schrodinger equation can be performed by the following 
choice of the wave functions ^f(p) as the eigenfunctions of J 2 , J3 and P: 

*(p) = -EV»i(py(n). (61) 
Pi 



Here the summa in r.h.s. of Eq. fl6"I|) runs over % — A, if P — (—)■', and over i — — , + if 
P = (— y +1 . Substituting this function into the new Schrodinger equation and collecting 
coefficients at bispinor harmonics (p A and (p° (or at (p~ and tp + ) one obtains the pair of 
coupled Rarita-Schwinger equations. In the matrix form they are: 

HV(p) = £ v|/(p), (62) 

where 

= MP)} (63) 

is two-component column wave function, 

1 f d 1 ,A 1 

and 

K= [K ik ] = [(i\L 2 + a 2 Z\k)] . (65) 

The form of 2x2 symmetric matrix K depends on both the parity and the tensor 
structure of interaction: 



K 



j(j + 1) + " 2 7 



2 v / i(i+ 1 ) 



(66) 
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for the parity P = (— ) J , and 



K 



(j+l)(j+2) + a 2 ( 7 - 



40+1) 



2j + l 



n r 1-S 2 3Vj(j+l) 
" 4(2j + l) + + l 



*»4(2j + l) + + l 



■> aj+i 



(67) 



for the parity P = (— ) J+1 . We note that in general case where ( ^ (including the cases 
of vector and gravitational interaction; c.f. Eqs.(|52|), (|33p) the operators £ and L 2 does 
not commute. Thus in calculating of K we chosen the ordering rule (0). The use of ([58]) 
or (|59| ) leads to off-diagonal elements of K which are singular at j = 1 (besides of j = 0). 
But there is no any physical reason for such a singularity. 

Using now an appropriate unitary (even orthogonal, to be sharp) transformation, the 
matrix K can be diagonalized, so that the coupled Rarita-Schwinger equations split into 
a pair of one- dimensional Schrodinger equations with effective Hamiltonians ^ of the 

form (p£2|) but with non-integer t(i, k y. 



(i,k) 



(i,k) 



(i,k), 



with 



K, 



(i,k) 



K u + K kk ± J(K U - K kk y + AKg 



(i^k), 



(68) 
(69) 



where i = A, k = if P = , and i = -, k = + if P = 

These equations are exactly solvable and lead to the mass spectrum (|T^)- ([l7|) but 
with another effective "principal quantum number": 



v 



v = P{n r ,j, P) =n r + i( itk ) + 1. 



(70) 



The calculation of effective "principal quantum number" v is straightforward by the use 
of Eqs. (|66|)-(f70|). Here we do not write down these rather cumbersome formulae. 



7 The confinement problem 

The Todorov recipe was observed on the systems with Coulomb-like interaction. Here 
we demonstrate that this rule appears useful for the construction of relativistic potential 
model of mesons. 

It is well known that the spectra of heavy quarkoniums are described satisfactory (and 
modulo spin effects) by means of the nonrelativistic potential model with the short-range 
Coulomb potential (0) and the long-range linear potential |T 



V y (r) = ar, (71) 

where a > is a constant. 

The description of light mesons needs the application of relativistic models. They 
frequently are built as single-particle wave equations with the vector short-range potential 
and the scalar long-range one |p2[ . Other models treat mesons as an extended objects or a 



composite two-quark relativistic systems. One of them which is concise and elegant, and 
which reflects principal features of the light meson spectroscopy, is the covariant oscillator 
model. Few versions of this model are given in Refs. 
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The nonrelativistic, single-particle and oscillator models appears related to one another 
by the Todorov recipe. Given the single-particle Klein-Gordon equation (|1|) with the scalar 
potential ((FTP and the vector one fl7|), this rule fixes unambiguously the form of the two- 
particle wave equation which is the relativization of the nonrelativistic potential 
model. 

If m a = and a = the wave equation reduces to the oscillator problem and yields 
the exact solution for the mass spectrum: 



M = 8a 



+ 2n r + | 



(72) 



The spectrum falls on the family of straight lines in the (M 2 ,£)-plane known in the hadron 
spectroscopy as the leading (for n r = 0) and daughter's (n r > 0) Regge trajectories. This 
structure and (£ + 2n r ) -degeneracy of the spectrum ([72] ) are characteristic of actual light 
meson spectra, if to neglect the rest mass contribution and fine spin effects. 

In the general case the equation ©-(0), (0), (|7T| ) is not exactly solvable. Here we use 



the oscillator approximation to estimate the spectrum for £ large. 

The substitution \l/(r) = ^(r)Y £ M (n) reduces the equation (|7lD to the form 

<(r) + [Q M - U Mi (r)]iP(r) = 0, (73) 

where 

U M i(r) = U M (r)+£(£+l)/r 2 . (74) 

The function Umi(t) nas a local minimum at some point r depending on M, £ and 
satisfying the condition 

U' m (r ) = 0. (75) 
Thus one can expand the potential (|74|) at the minimum, 

U M i{r) » U Mi (r ) + ]fl u (r )(r - r ) 2 , (76) 
and search a solution of this oscillator problem. A quantization condition then reads: 

Qm - U M£ {r ) = ^\U'l u {r Q ){2n T + 1). (77) 

If n r ~ 1, the approximate solution differs exponentially little from the exact solutions of 
the problem. 

Eqs. (|75|) and (|77|) form the set of algebraic equations with r and M to be found. 
Solving this set by a power series in £ leads to the asymptotic formulae: 



and 



M 2 



+ 2n r + l- V2a\ + 2 (m\ + m 2 + \f2m x m<^ + 0{r l ). (79) 

The latter represents the spectrum of the system at £ ^> n r (but provides a good approxi- 
mation even if £ ~ 2^-3). As to compare this formula to the Eq. (0), the influence of the 
short-range interaction and non-zero rest masses result in the parallel shift of the family 
of Regge trajectories as a whole. Considering the constants mi, m 2 (and, possibly, a) as 
adjustable parameters one can obtain trajectories for different families of light mesons. 
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Up to now we neglected a mass splitting due to a spin interaction. The majority of 
attempts to describe spin effects in hadron spectroscopy concerns with the heavy quark 
systems which can be treated as weakly relativistic systems. As usual, one takes the 
generalized Pauli-Fermi Hamiltonian with long-range scalar potential and short-range 
vector one (the potentials (fZlj) , (|7|) in our case) and corresponding spin corrections (|I9|), 



(|20D treated perturbatively [17|, |24| . But this scheme can fail when considering light 
mesons as corresponding to strongly relativistic domain M 3> mim- First of all we note 
that, as it follows from Eqs. (]7q) , (|79|) , in this domain r 2 ~ £ and M 2 ~ £. Thus r ~ M, 
i.e., the radius of meson is proportional to its mass; here we took into account that the 
wave function at £ large is localized around tq. Then we have rough estimates: 



V s = ar~M, W s ~£/r~M V s ~ W s , 

V v = -a/r ~ M~\ W v ~ £/r 3 ~ M" 1 K ~ W v . (80) 

The spin corrections appears to be of the same order as the nonrelativistic potentials. 



But actual spin effects in light meson spectra are small. Moreover, the operators (19) 
and (P0|) was deduced within the perturbation theory |TB|, |T7|, so they should satisfy 
inequality W <C V by construction. Second, the operators fll9"|) and ( p0|) are divergent 
if mi or/and m,2 — > 0. Consequently, divergent terms appear in the mass spectrum (in 
contrast to the case of Coulomb-like interaction where the spectrum is not singular if 
mi or/and ni2 vanish). 

The possible way to avoid these two problems is the following modification of operators 
(|19"D and (|(J): the rest masses of quarks involved in these operators should be replaced 
by the "constituent" masses: 

M 2 + m 2 — mi 



m a ^ M a = ^m 2 a + Q M = ^ * a = 1,2, a = 3-a (81) 

possessing the properties 1) .\/, • M> M, 2) M a w M/2 if M > m^), and 3) M a -> m a 
in the nonrelativistic limit M — > m + (we suppose that M 2 > |m 2 — m 2 ,!). The modified 
operators W are regular at mi or/and rri2 — > 0, they became small as ~ ^(m + /M) 2 in 
strongly relativistic domain M ^> mu2), and reduce to W in weakly relativistic domain 
M w m + . 

Now let us estimate the effect of spin corrections W s to the scalar potential (f7I|) 
(the effect of vector interaction is minor, as it follows from the estimates in ( |S0"D and 
the paragraph above). For this purpose we modify the scalar potential V s = V s + W s 
and substitute it (instead of V s ) into the effective potential Um- Then we diagonalize the 
modified effective potential U m{t) and apply the oscillator approximation to the resulting 
pair of split quasipotential equations. Finally, we come to the asymptotically linear Regge 
trajectories ( |79] ) with the same slope parameter 8a, but each trajectory splits into three 
ones by parallel shift 

o [ 0, l = j, 
AM 2 = { (82) 
[ ±2a, £ = j±l. 

Note that the mass splitting does not depend on £. Qualitatively this result correlates 
with actual light meson spectra as well as with theoretical results following from the string 



models |25|, [26| . 

As in the spinless case, the model is exactly solvable if m a = and a = 0. The 
mass spectrum is determined by M 2 = 4ax, where x is the positive solution of the 
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transcendental equation: 



x 2 - 2{2n r + 1)x-k = y / (2£+l) 2 x 2 + k 2 . (83) 
which, in turns, reduces to a cubic algebraic equation, and 

'o, e = j, s = o, 

1, i = J, s = 1, 

K=l (84) 

e + i, £ = j + i, s = i, 

t -£, I = j - 1, 8 = I. 

Here we do not write down an explicit form of the mass spectrum. 

Of course, the minimal set of adjustable parameters makes the present potential model 
too poor to provide a sharp fit to an experimental data. In particular, the mass splitting 
Q82D is unambiguously fixed by the typical hadron scale 8a. This value (0% and 25% of 
8a) contradicts to some estimates of actual experimental data (about 5-^6%; see [f2~6|j). 

In another version of the model we include operators W similarly to the case of 
Coulomb-like interaction (see Eq. (P2"|)), i.e., by the following modification of the effective 
potential: 

U M {r) -> U M {r) = U M {r) + 2m M W s + 2E M W V . (85) 

In the strongly relativistic domain the present potential model leads to the same asymp- 
totic spectrum (|79|) as the spinless model does. The mass splitting is small. It can be 
calculated perturbatively and fitted to experimental data by adjusting the parameters m a 
and a. 



8 Summary 

In the present paper we have constructed the family of simple quantum-mechanical models 
which describe two-fermion relativistic systems with different Coulomb-like and confining 
interactions. We have embodied in these models some theoretical experience of studies 
in the relativistic two-body problem. We started with the quasipotential equations de- 
scribing the scalar and vector interaction of two spinless particles. As it was observed by 
Todorov m [/J, these equations are solvable, and they have simple single-particle structure. 
Similar features are characteristic of the family of RDIT models W\ 13[] which describe 
an arbitrary relativistic Coulomb-like interaction including the gravitation. We generalize 
these equations to the case of two-fermionic systems. For this purpose we use the oper- 
ators known from QFT which describe spin-dependent corrections to the scalar, vector 
and gravitational interacting ||, [TB], [T7|, [T3| . We first treat these operators perturbatively 
and obtain the spectrum of muonium and its scalar and gravitational counterparts with 
accuracy up to a 4 . Then we modify the equations in such a way that they become exactly 
solvable and yield correct (within the same accuracy) mass spectra. These equations can 
be useful for accounting higher than a 4 (say, radiative) corrections to a particle interaction 
by means of the first order perturbation theory. 

Also we have demonstrated that the Todorov recipe of constructing two-body equa- 
tions permits a straightforward application to the case of confining interaction. The 
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generalization to spinning particles has been proposed too. In the weakly relativistic do- 



main this equation reduces to well-known potential model [I7[] which is appropriate to the 
description of heavy mesons. In the strongly relativistic limit it yields a mass spectrum 
which reproduces qualitatively light-meson experimental data. 
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